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1.  Foreword 


Relatively  little  is  known  about  helicoptei  stability  during  high-load-factor  maneuvers  under 
unsteady  conditions  (Refs  1-3)  These  maneuvers  represent  short-duration  dynamics  that 
lasts  long  enough  to  stall  the  rotor  but  not  long  enough  to  reach  a  steady  state  They  also 
represent  operations  at  the  extremes  of  the  flight  envelope  wher  e  the  peak  loading  goes  well 
beyond  the  steady  thrust  limit,  dynamic  stall  occurs  and  often  a  design  condition  is  reached 
(Refs  2-4)  (Routine  transient  maneuvers,  which  do  not  last  long  enough  to  stall  the  rotor  as 
well  as  steady  turns,  which  can  be  treated  on  the  basis  of  Floquet  theory,  are  not  addressed 
here ) 

Presently,  helicopter  stability  is  computed  from  state-space  models  and  from  numerically  or 
experimentally  generated  responses  by  a  host  of  algorithms,  ranging  from  the 
moving-block  analysis  to  sparse-time-domain  technique  to  generalized  fast-Floquet  theory 
(Refs  5,  6).  No  matter  how  these  algorithms  are  designated,  it  is  Floquet  theory  that 
under  lies  all  of  them  (Refs  5,  6)  Floquet  theory  is  a  linear  stability  theory,  and  it  provides 
asymptotic  stability  under  periodic-response  conditions  Aggressive  maneuvers  of  short 
duration  do  not  reach  a  steady  state  What  about  stability  on  the  basis  of  Lyapunov 
exponents,  which,  after  all,  is  a  generalization  of  linear  stability  based  on  Floquet  theory? 
For  example,  the  Lyapunov  exponent  reduces  to  the  real  part  of  the  Floquet  exponent  for  a 
periodic  orbit  and  provides  a  means  of  quantifying  stability  of  steady-state  response  that  need 
not  necessarily  be  periodic  (e  g  quasi-per iodic  orbit,  Refs  1,  2,  7)  The  Lyapunov  exponent, 
however,  requires  averaging  of  response  data  over  a  long  time  Simply  put,  neither  the 
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Lyapunov  exponent  approach  nor  the  Floquet  approach  a  fortiori  provides  a  means  of 
tr  eating  stability  of  aggressive  maneuver  s  of  short  dur  ation 

The  task  of  treating  stability  of  aggressive  maneuvers  of  short  duration  moves  into  an 
uncharted  territory  of  nonlinear  stability  under  unsteady  conditions  The  STIR  work  does  just 
that  Specifically,  it  explores  how  the  emerging  finite-time  Lyapunov  exponents  (FILE) 
provide  the  key  to  determine  stability  and  thus  to  provide  a  framework  to  calculate  stability 
(Refs..  2,  7-9)  This  work  does  not  address  chaos;  specifically,  it  addresses  how  the  largest 
FTLE  can  be  exploited  to  quantify  the  stability  margin  of  the  least  stable  mode  during  such 
maneuvers 

To  help  appreciate  the  relevance  of  the  FTLE  approach,  this  report  begins  with  a  passing 
reference  to  the  RAFI-66  Comanche  helicopter  lesson  and  to  maneuvers  under  steady  and 
unsteady  conditions  Then  it  comes  to  FTLE  as  a  stability  measure  for  aggressive  maneuvers 
of  short-duration  After  presenting  conclusions,  it  also  mentions  the  continuing  study  of  the 
FTLE  approach,  including  its  application  to  a  comprehensive  set  of  numerically  and 
experimentally  generated  data;  such  a  continuing  study  should  bring  out  the  strengths  and 
weaknesses  of  this  appr  oach 


2.  Aggr  essive  Maneuvers 

The  US  Army  has  been  emphasizing  the  development  of  air  vehicles  with  increasingly 
stringent  maneuvering  requirements  for  their  combat  missions  such  as  scout-attack  (Ref!  4) 
These  maneuvers  are  typically  aggressive  maneuvers,  mostly  of  short  duration,  and  often 
they  represent  a  design  condition 
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Fig  1  Helicopter  Schematic  in  a  Diving  Bank-Angle  Turn 

Of  particular  significance  is  the  lesson  from  the  RAH-66  Comanche  helicopter  development 
This  vehicle  was  found  to  reach  a  steady  condition  during  aggressive  diving  turns  (Fig.  1); 
also,  see  next  section  for  details  on  steady  and  unsteady  conditions  Although  this  STIR  work 
addresses  unsteady  conditions,  the  stability  issue  that  plagued  the  Comanche  development 
merits  some  elaboration 

Therein,  inadequate  stability  margin  or  blade  damping  was  discovered  during  flight  tests  of 
diving  turns  Specifically,  the  vehicle  was  found  to  be  vulnerable  to  air  resonance  instability 
in  which  the  blade  regressing  lag  mode  (RLM)  frequency  coalesced  with  the  fuselage  roll 
mode  frequency,  and  in  turn  the  blade  inplane  damping  sharply  declined  with  increasing 
thrust  level  and  was  later  augmented  by  an  active  feedback  controller  (Ref  10) 
The  state-of-the-art  prediction  codes  failed  to  predict  the  decline;  in  fact,  conventional 
wisdom  has  it  that  stability  is  not  an  issue  during  maneuvers..  This  failure  demonstrated  that 
stability  during  aggressive  maneuvers  is  a  crucially  important  and  yet  umesolved  issue  in  the 
Army’s  development  of  air  vehicles  for  combat  missions 
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3,  Steady  and  Unsteady  Conditions 

During  forward  flight,  it  is  not  difficult  to  distinguish  between  steady  and  unsteady 
conditions  and  in  turn  to  conceptualize  perturbation  from  a  periodic  orbit  as  required  by 
Floquet  theory  By  contrast,  “the  distinction  between  the  steady  and  unsteady  maneuvers  is 


sometimes  difficult  to  make”  (Ref  4);  after 
accelerated  flight 


all,  maneuvers  do  involve  some  form  of 
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Fig  2  “Load  factor  time  histories  in  two  level  bank-angle  turns,  illustrating  maneuvers  steadiness 
Duration  time,  Td,  indicates  portion  of  time  history  where  load  factor  is  within  2%  of  maximum 
value  a)  Counter  8539;  b)  Counter  8826  ” 


Fig  3  “Load  factor  time  histories  in  two  diving  bank-angle  turns,  illustrating  maneuvers  steadiness 
Duration  time,  Td,  indicates  portion  of  time  history  where  load  factor  is  within  2%  of  maximum 
value  a)  Counter  11668;  b)  Counter  11683  ” 
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For  a  better  appreciation  of  this  subtle  distinction,  it  is  instructive  to  study  Figs  2  and  3, 
which  show  load-factor  time  histories  for  the  level  bank-angle  and  diving  bank-angle  turns 
from  the  UH-60A  airloads  program  (Ref  4)  As  seen  from  Fig  2a  (counter  8539),  the 
maximum  load  factor  nz  is  1  52g  and  the  horizontal  dotted  line  is  98%  of  the  maximum  or 
1  49g  The  maneuver  point  duration  Td  is  defined  for  the  period  during  which  the  load  factor 
nz  exceeds  the  98%  level,  that  is,  within  2%  of  maximum  value;  this  duration  is  taken  as  a 
measure  of  the  steadiness  of  the  maneuver  (Ref  4)  While  the  maneuver  point  duration 
Td  =  11  seconds  in  Fig  2a,  Td  =  4  2  seconds  in  Fig  2b,  about  4  times  longer  Qualitatively, 
the  same  type  of  variation  is  depicted  by  the  diving-turn  case  in  Fig  3a  (counter  1 1668)  and 
in  Fig  3b  (counter  1 1683)..  Although  the  identification  of  Td  and  the  corresponding  measure 
of  maneuver  steadiness  are  arbitrary,  Figs  2  and  3  do  show  that  these  maneuvers  basically 
fall  into  two  categories:  In  the  first,  the  maneuvers  do  exhibit  steady  conditions;  that  is,  the 
maneuver  point  duration  Td  contains  multiple  consecutive  rotor  revolutions,  which  can  be 
considered  a  steady  (periodic)  condition  from  Floquet  theory  stand  point  In  the  second,  such 
a  steady  condition  is  not  exhibited,  and  the  Floquet  theory  is  no  longer  applicable  or  its 
applicability  is  nebulous  at  best 


4.  Backgr  ound  to  FTLE 

Ihe  Lyapunov  exponent  is  the  most  fundamental  measure  or  metric  of  stability,  Its 
computation  requires  time  averaging  of  the  response  vector  over  a  long  time  (Refs  7-9) 
In  particular,  it  reduces  to  the  real  part  of  the  Floquet  exponent  for  the  special  case  of  a 
periodic  orbit  The  theoretical  basis,  computation  and  application  of  FTLE  represent  an 
emerging  and  active  area  of  research  (Refs  9,  1 1,  12)  FTLE  is  the  finite-time  version  of  the 
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Lyapunov  exponent,  and  it  provides  a  rigorous  means  of  quantifying  the  finite-time  or  local 
stability  for  that  duration  of  the  orbit  and  for  that  initial  state  In  the  next  section,  stability 
formulation  based  on  FTLE  is  presented  with  a  passing  reference  to  the  Lyapunov  exponent, 
in  particular,  to  the  connection  between  F  TLE  and  the  Lyapunov  exponent  As  a  background 
to  FTLE,  however,  a  brief  account  of  the  widely  used  generalized  Floquet  approach  is 
included  in  the  sequel  (Ref  5)  Like  the  generalized  Floquet  approach,  the  FTLE  formulation 
is  not  restricted  to  state-space  modeling 

4.1  Gener  alized  Floquet  Appr  oach 

Classical  Floquet  theory  applies  to  state-space  models  For  a  model,  say  with  M  states,  it  is 
required  to  perturb  or  excite  each  of  the  M  states  about  a  periodic  response,  one  state  at  a 
time,  to  generate  the  M  x  M  perturbation  matrix  P  at  some  starting  time,  say  at  t  =  to-  It  also 
requires  that  each  of  the  corresponding  M  responses  are  measured  at  the  end  of  one  period, 
say  at  t  =  to  +  2tc,  to  generate  the  M  x  M  response  matrix  R  The  Floquet  transition  matrix 
(FTM)  connects  P  and  R: 

O^Imxm  =[ftm]  MxM  MmxM  _ (1) 

The  generalized  Floquet  theory  provides  an  analytical  basis  to  predict  the  stability  of  models 
without  a  state-space  representation  and  with  numerically  or  experimentally  generated 
response  histories  Algorithmically,  it  also  connects  P  and  R,  but  it  permits  perturbation  of 
an  arbitrary  number  of  states  and  generalizes  the  above  equation  (Refs  1,  2,  5): 

M i x j  =[FTM]IxI  [P]IXJ  _ (2) 

The  number  of  rows  I  and  the  number  of  columns  J  depend  upon  the  number  of  states 
perturbed,  the  number  of  additional  periods  and  time  delays  used  in  generating  the 
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time-shifted  or  pseudo  states  such  as  Xj  (t  +  At)and  xs  (t  +  2At),  where  X;(t)  is  a  state  and 
At  is  the  time  delay  or  time  shift  Now,  the  FTM  is  generated  by  the  singular  value 
decomposition  of  P  and  R,  and  then  by  P*,  the  generalized  inverse  of  P  That  is,  the  FTM  is 
approximated  in  the  least-squares  sense  Stated  explicitly, 

[ftm]=[r]  [p]  *  _ <3) 

5.  FTLE  Approach 

While  the  Floquet  exponents  describes  the  stability  of  a  periodic  response,  the  infinite-time 
or  global  Lyapunov  exponent  is  a  generalization  of  the  Floquet  exponent  in  that  it  describes 
the  stability  of  an  orbit  that  is  not  constrained  to  be  periodic  Similarly,  FTLE  describes 
stability  of  an  orbit  that  does  not  last  long  enough  to  reach  a  steady  state,  and  thus  it  provides 
the  basis  for  the  framework  to  calculate  the  stability  This  framework  involves  three 
sequential  computational  blocks:  i)  constructing  a  pseudo-state  space;  ii)  generating  a  series 
of  partial  derivative  or  Jacobian  matrices  and  in  turn  forming  the  positive  symmetric 
Oseledec  matrix,  and  iii)  diagonalizing  this  Oseledec  matrix  by  the  recursive 
QR  decomposition  (Ref  9) 

5.1  Pseudo-state  Space 

The  pseudo-state  space  construction  begins  with  N-point  time  series  of  the  response: 

{xi,  x2,  x3 . Xj...  xN},  where  Xj  is  a  scalar  measurement  at  discrete  time  t  =  tj.  Now  it  is 

required  to  introduce  four  important  parameters: 

At  =  data  sampling  interval  (assumed  constant)  =  f+i- 1\ 
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m  =  embedding  dimension 
X  -  time  lag  =  (integei)  At 

J  =  delay  expressed  in  number  of  data  points  —  x  /  At 

With  these  ingredients,  an  Mxm  matrix  [X]  with  M  vectors  is  generated: 

X,  =  [x„x1+J, . ,  xi+(m_i)j] 

Xj  =  [X  j>X2J  5  j  XmJ  ] 

XM  =  [XM>XM  +  J> ■"»  XM  +  (m-l)j] 

where  X*  is  the  reconstructed  state  of  the  system  at  t  =  tj 

It  is  verified  that  N  =  M  +  (m-1)  J  The  local  Lyapunov  exponents  are  computed  in  the 
embedded  space  of  these  M  vectors 

5.2  Towar  d  Computing  FTLE 

Say  m0  Lyapunov  exponents  are  computed.  With  k  representing  the  starting  time  step 
t  =  h  =  k,  the  stipulation  is  that  the  response  vector  X(k)  is  governed  by  the  mapping 
(Refs.  2,  8,  9): 

X(k)  — » X(k  +  l)  =  F(x(k)) 

Roughly  speaking,  stability  shows  how  perturbation  AX(k)  to  X(k)  grows  or  decays  with 

time  Accordingly,  to  connect  A(k)  with  A(k+1),  set 

X(k+l)+A(k  +  l)  =  F(X(k)+A(k)) 

=  F(X(k))  +  DF(X(k))A(k) 

In  other  words, 

A(k  +  l)  =  DF(X(k))  A(k) 
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where  DP  (X(k))  represents  the  m0x  m0  Jacobian  or  partial  der  ivative  matrix  such  that 

DF(X) .  =-5^;  i  =  1,  2,  ..  ,  mo,  and  j  -  1,2,  , 
dX. 

J  m0 

It  is  instructive  to  observe  that  whenX(k  +  2rc)  =  X(k) ,  the  stability  is  governed  by  Floquet 

theory  After  L  time  steps,  it  is  verified  that 

A(k  +  L)  =  DF(X(k  +  L-l))  DF(X(k  +  L-2))  DF(X(k))  A(k) 

=  DF[  (X(k))  A(k)  . (4) 

By  definition,  the  positive  symmetric  Oseledec  matrix  after  L  time  steps  is  given  by 

os(x,l)  =  [df1  (x)]  [df1(x)]1  - (5) 

The  matrix  has  real  eigenvalues;  their  logarithms  are  called  the  finite-time  Lyapunov 
exponents  and  are  represented  by  A*  (x,  L),  X\  (x,  L)  being  the  largest  Although  the  Oseledec 
matrix  is  ill-conditioned,  the  recursive  QR  decomposition  provides  a  means  of  computing  the 
eigenvalues  It  is  instructive  to  observe  that 

lim^i(x,L)=Xi,  \>X2>  >Xm  - (6) 

L— 

where  the  logarithm  of  A.;  is  the  conventional  Lyapunov  exponent  For  regular  motion,  X\  <  0 
and  for  chaos,  Ai  >  0  Hereafter,  the  treatment  is  restricted  to  stable  regular  motion  and  to  the 
largest  FTLE,  which  is  simply  designated  as  Ail  (<  0) 

5.3  FTLE  as  a  Stability  Measur  e 

Stability  means  the  exponential  rate  of  growth  or  decay  with  time  of  two  neighboring 
trajectories;  this  is  shown  schematically  in  Fig  4  Therein,  5  (t)  represents  the  separation 
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between  the  trajectory  x  (t)  with  initial  state  xG  and  the  peituibed  trajectory  x  (t)  +  5  (t)  with 
initial  state  x0  +  S0  After  L  time  steps,  the  separation  5  (t)  can  be  r  epr  esented  as 

|8(t)|HKI|el'“  - (7) 

With  L— >  oo,  Xil— >  Ai,  and  for  a  periodic  trajectory  with  period  T, 

^i=YInN  - <s) 

where  pi  is  the  largest  characteristic  multiplier  or  the  largest  eigenvalue  of  the  FTM,  given 
by  Bq  (3)  For  a  stable  system,  it  is  instructive  to  observe  that  -A-i  (>  0)  provides  the  stability 
margin  of  the  least  stable  mode 


Fig  4  Schematic  of  Separation  between  two  Neighboring  Trajectories 

For  rotorcraft  systems  designed  for  aggressive  maneuvers,  the  stability  issue  often  centers  on 
how  stable  is  the  least  stable  mode,  that  is,  on  the  adequacy  of  the  stability  margin  from 
handling  qualities  perspective  Thus,  what  is  required  is  a  means  of  quantifying  the  stability 
margin..  Recall  that  -A*  (>  0)  quantifies  the  stability  margin  for  the  least  stable  mode;  it  is 
evaluated  from  Floquet  theory  for  a  periodic  trajectory  and  from  Lyapunov-exponent  method 
for  a  trajectory  that  is  not  constrained  to  be  periodic  Furthermore,  AIL— ►  h  with  L— >  co; 
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that  is,  Xu  selves  as  the  finite-time  version  of  X{  All  this  points  to  the  feasibility  of 
using  -X\i  (>  0)  as  a  sensible  means  of  quantifying  the  stability  margin  of  the  least  stable 
mode  for  aggr  essive  maneuvers  of  short  duration 

6.  Conclusions 

The  investigation  toward  developing  a  framework  to  calculate  helicopter  stability  during 
aggressive,  short-duration  maneuvers  has  been  presented  This  investigation  exploits  the 
largest  finite  time  Lyapunov  exponent  as  a  measure  of  stability  margin,  and  it  represents  the 
first  attempt  toward  such  a  development  from  first  principles  It  is  still  in  an  exploratory 
stage  The  next  section  gives  an  account  of  its  continuation  under  a  three-year  ARO  grant 
(May  1,2007 -April  30,  2010) 


7..  Continuing  Study 

Two  major  objectives  of  the  continuing  study  ate  1)  to  further  refine  the  analytical  basis  of 
the  developed  framework  to  calculate  the  stability  and  2)  to  investigate  the  computational 
issues 

The  computation  of  the  largest  finite  time  Lyapunov  exponent  XlL  represents  an  active 
research  area  and  is  computationally  demanding.  One  reason  is  that  the  Oseledec  matrix 
given  by  Eq  (5)  is  ill-conditioned  Accordingly,  the  recursive  QR  decomposition  will  be 
investigated  to  reliably  compute  LiL  (Ref  9) 
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Finally,  the  database  from  the  UH-60  aii loads  program  will  be  used  to  identify  aggressive 
maneuvers  of  short  duration  and  then  to  calculate  the  stability  margins  for  these  maneuvers 
This  should  bring  out  the  strengths  and  weaknesses  of  the  framework  to  calculate  the  stability 
of  aggressive,  short-duration  maneuvers  on  the  basis  of  the  largest  finite-time  Lyapunov 
exponent 
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